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(57) A method and system for thermal-dynamic 
modeling and performance evaluation of a nuclear Boil- 
ing Water Reactor (BWR) core design is presented. A 
data processing system is used to execute specific pro- 
gram routines that simultaneously simulate the thermal 
operating characteristics of fuel rods within the reactor 
during a transient operational condition. The method 
employs a multi-dimensional approach for the simula- 
tion of postulated operational events or an anticipated 
operational occurrence (AOO) which produces a tran- 
sient condition in the reactor — such as might be caused 
by single operator error or equipment malfunction. 



Based on a generic transient bias and uncertainty in the 
change in critical power ratio (ACPR/ICPR), histograms 
of fuel rod critical power ratio (CPR) are generated. Ul- 
timately, the operating limit minimum critical power ratio 
(OLMCPR) of the reactor is evaluated from a histogram 
of probability calculations representing the number of fu- 
el rods subject to a boiling transition (N RSBT) during the 
transient condition. The histogram may be readily dis- 
played by the data processing system and used to sta- 
tistically demonstrate an OLMCPR compliance of the re- 
actor core design with USNRC regulations. 
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Description 

[0001] The present invention relates generally to methods for evaluating nuclear power core operation, and more 
particularly to an improved method and apparatus for determining an operating limit minimum critical power ratio (OL- 
MCPR) so as to effectuate increased efficiency and operation of Boiling Water Reactors (BWR). 
[0002] A Boiling Water Nuclear Reactor (BWR) generates power from a controlled nuclear fission reaction. As shown 
in Figure 1 , a simplified BWR includes a reactor chamber 101 that contains a nuclear fuel core and water. Generated 
steam is transferred through pipe 1 02 to turbine 1 03, where electric power is generated, then water returns to the core 
through pipe 1 04. As shown in Figure 2, the core 201 is made of approximately five hundred (500) bundles 202 of fuel 
rods arranged in a specific manner within the reactor core. As shown in Figure 3, each bundle 301 contains roughly 
one hundred (100) fuel rods 302. Water in the core surrounds the rods. Heat generated by a nuclear reaction is trans- 
ferred from the rods to the water circulating through the core, boiling some of the water. The heat generated in the core 
is carefully controlled to maintain safe and efficient reactor operations. 

[0003] In a Boiling Water nuclear Reactor there are basically three modes of heat transfer that must be considered 
in defining thermal limits for the reactor: (i) Nucleate boiling, (ii) transition boiling and (iii) film boiling. "Nucleate boiling- 
is the preferred efficient mode of heat transfer in which the BWR is designed to operate. "Transition boiling" is manifested 
by an unstable fuel rod cladding surface temperature which rises suddenly as steam blanketing of the heat transfer 
surface on the rod occurs, then drops to the nucleate boiling temperature as the steam blanket is swept away by the 
coolant flow, and then rises again. At still higher fuel rod/bundle operating powers, "film boiling" occurs which results 
in higher fuel rod cladding temperatures. The cladding temperature in film boiling, and possibly the temperature peaks 
in transition boiling, may reach values which could cause weakening of the rod cladding and accelerated corrosion. 
Fuel rod overheating is conservatively defined as the onset of the transition from nucleate boiling to film boiling. The 
conventional basis for reactor core and fuel rod design is defined such that some "margin," accommodating various 
design and operational "uncertainties," is maintained between the most limiting operating condition and the transition 
boiling condition at all times for the life of the core. 

[0004] The onset of transition boiling can be predicted by a correlation to the steam quality at which boiling transition 
occurs— called the "critical quality." Steam quality can be readily measured and is generally a function of measuring 
distance above the boiling boundary (boiling length) for any given mass flow rate, power level, pressure and bundle 
flow geometry among other factors. A "critical power" is defined as that bundle power which would produce the critical 
quality of steam. Accordingly, a "critical power ratio" (CPR) is defined as the ratio of the critical power to the bundle 
operating power at the reactor condition of interest. CPR is descriptive of the relationship between normal operating 
conditions and conditions which produce a boiling transition. CPR is conventionally used to rate reactor design and 
operation. To assure a safe and efficient operation of the reactor, the CPR is kept above a prescribed value for all of 
fuel assemblies in the core. Reactor operating limits are conventionally defined in terms of the most limiting fuel bundle 
assembly in the core— defined as the "minimum critical power ratio" (MCPR). Reactor operating limits are often stated 
in terms of MCPR. 

[0005] In nuclear power generation engineering, it is widely recognized that there is a possibility, however small, that 
the occurrence of a reactortransient event combined with the various "uncertainties" and tolerances inherent in reactor 
design and operation may cause transition boiling to exist locally at a fuel rod for some period of time. Accordingly 
MCPR operating limits are conventionally set in accordance with the United States Nuclear Regulatory Commission 
(USN RC) design basis requirement that transients caused by a single operator error or a single equipment malfunction 
shall be limited such that, taking into consideration uncertainties in the core operating state, more than 99.9% of the 
fuel rods are expected to avoid boiling transition during that error or malfunction. A safety limit minimum critical power 
ratio (SLMCPR) is defined under current USNRC requirements as the MCPR where no more than 0.1% of the fuel 
rods are subject to boiling transition. The corresponding OLMCPR describes the core operating conditions such that 
the MCPR is not lower than the SLMCPR to a certain statistical confidence. 

i. An Ideal Approach 

[0006] In principle, the OLMCPR could be calculated directly such that for the limiting anticipated operational occur- 
rence (AOO), less than 0.1% of the rods in the core would be expected to experience boiling transition. This approach 
is described in U.S. Patent No. 5,912,933, by Shaug et al. The process involved is shown in Figure 4, which depicts 
histogram 400 of rod CPR values 401 verses number of rods 402 at the specific CPR value. While the CPR value is 
usually associated with a fuel bundle, it actually refers to the limiting rod in a bundle. Each rod in the bundle has a CPR 
value that is determined by the local power distribution (TIP uncertainty) and relative position of the rod within the 
bundle (R-factor). The lowest CPR of any one rod in the bundle is used to characterize the CPR for the entire bundle. 
[0007] The CPR 401 for a given rod has an associated probability distribution function (PDF) which reflects the 
uncertainties in its determination. The PDF may be determined experimentally and is shown as an Experimental Critical 
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Power Ratio (ECPR) distribution 410. Thus, if a nominal CPR value (411) is 1.0, then the PDF 410 of probable actual 
CPR values range from 0.90 to 1.10. The variability in the rod CPR values is due to uncertainties in the initial rod 
condition, i.e., uncertainties in the measurements of parameters at the reactor operating state (core power) and in the 
modeling of derived parameters (power distribution). 

5 [0008] To take the effect of a transient event on the CPR values into account, a safety margin is introduced to CPR 
values by shifting the acceptable nominal CPR value 405 for the lowest rod CPR to a larger CPR value, i.e., 1 .25. The 
ECPR histogram distribution 403 for the lowest CPR rod is thus shifted such that the entire CPR histogram is above 
a CPR value of 1 .20, and well above a CPR value of 1 .0. Moreover, the nominal CPR values 407 for rods other than, 
the lowest CPR rod are above the nominal CPR value, e.g., 1 .25, of the lowest CPR rod. 

w [0009] During a transient in rod operation, the histogram 407 of rod CPRs shifts to the left to lower CPR values, 
resulting in the histogram 408. With this shift, the "nominal" CPR value 406 during the transient is at the point, e.g., 
1 .05, where the minimum CPR value is reached during the transient. The limiting rod will have an associated PDF 404, 
which includes both the uncertainties in the initial rod conditions and "transient uncertainties." The maximum change 
in critical power ratio during the transient (ACPR 409) includes uncertainties in the modeling of the transient, both the 

75 physical models and plant parameters. Ideally, this transient ACPR 409 and associated OLMCPR would be generated 
as shown in Figure 5, and described as follows: 

Step 1 : Assume a set of base core operating conditions using the parameters to run the plant generates a core 
MCPR equal to the OLMCPR, as shown by block 501 . 

20 

Step 2: Using the parameters, such as core power, core flow, core pressure and others, that predict a general 
bundle CPR set forth in block 506, determine the MCPR for each bundle in the core, as shown by block 502. 



Step 3: Using parameters, such as rod placement within each bundle and rod power, that change each bundle 
25 CPR into individual rod CPR values set forth in block 507, determine the MCPR for each rod in the core, as shown 

by block 503. 



Step 4: Using the ECPR probability distribution, generated by Equations 1 and 2, set forth in block 508, determine 
the percentage of NRSBT in the core by summing the probabilities of each rod in the core that is subject to boiling 
30 transition, as shown by block 504. This summation is shown by Equation 3. 



ECPR = (CriticatPower Predicted by Correlation) (Equation 1 ) 

(Measured CnticalPower) 



35 



00 '-u 2 



40 

(Equation 2) 

45 NRSTB (%) = x ]T [P { x (1 Rod)] 



(Equation 3) 

Step 5: Vary the parameters set forth in blocks 506 and 507 for a set number of Monte Carlo statistical trials, as 
shown by block 505. Compile the statistics from all the trials from steps 2 through 4 to generate a probability 
distribution of NRSBT. 



55 step 6: Compare the value of NRSBT percentage to 0.1%, as shown in block 509. If the percentage is greater 

than 0.1 %, reset the core parameters to different initial conditions in order to comply with the USNRC regulations, 
as shown in block 510. Similar to Step 1 and block 501, the new initial conditions are assumed to generate an 
OLMCPR. The determination of NRSBT restarts and loops until the value of NRSBT is equal to 0.1%. Similarly, if 
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the percentage is less than 0.1%, the core parameters are reset to increase the value of N RSBT in order to operate 
the core more efficiently or with fewer effluents. 

Step 7: If the percentage of NRSBT equals 0.1%, the assumed value of OLMCPR, which equals core MCPR, 
5 complies with the USNRC regulations, as shown by block 511 . Accordingly, the operating core conditions are set 

as the assumed parameters. 

[001 0] Because of the computational difficulties and the need to evolve efficient algorithms, the ideal process outlined 
above has not been followed. The currently approved process and the new approach to determining the OLMCPR are 
10 described in the following sections. 

ii. The USNRC Approved Approach 

[0011] In the current process, the OLMCPR determination is divided into primarily two steps, as shown by Figure 6. 

is Using a process similar to the ideal process, first the SLMCPR is determined so that less than 0.1% of the rods in the 
core will be expected to experience boiling transition at this value. In other words, 99.9% of the fuel rods in the core 
will be expected to avoid boiling transition if the MCPR in the core is greater than SLMCPR. Second, the OLMCPR is 
then established by summing the maximum change in MCPR (ACPR 95/95 ) expected from the most limiting transient 
event and the SLMCPR. Since Figure 6 is similar to the Figure 4, only a brief description of its components follows. 

20 Histogram 600 shows the number of rods at a specific CPR value 602 versus the corresponding CPR value 601 . The 
histogram 608 results with the lowest CPR rod 607 at a value of, e.g., 1 .05, which equals the SLMCPR 603. Limiting 
rod distribution 606 shows the uncertainty in determination of the limiting CPR rod 607. Similar to the ideal process, 
the SLMCPR 603 is determined when the percentage of NRSBT is equal to 0.1%. However, unlike the ideal process, 
the current process is unable to fully predict and measure certain parameters, such as the power distribution within 

25 each bundle and the power distribution along each rod. Thus, the uncertainties in calculating the SLMCPR do not allow 
equating the OLMCPR with the SLMCPR. Accordingly, an error factor, ACPR 95/95 605, is linearly added to the SLMCPR 
603 to determine the OLMCPR 609. ACPR^s 605 conservatively corrects for the inherent limitations in the calculation 
of the SLMCPR 603. 

[0012] Using the currently approved process, the OLMCPR 609 is generated as shown in Figure 7, and described 
30 as follows: 

Step 1 : Assume a set of base core operating conditions using the parameters to run the plant generates a core 
MCPR equal to the SLMCPR as shown by block 701 . 

35 Step 2: Using the parameters, such as core power core flow, core pressure, bundle power and others, that predict 

a general bundle CPR set forth in block 706, determine the MCPR for each bundle in the core as shown by block 
702. This process step has large uncertainties in predicting the bundle power, biasing the calculations. 

Step 3: Using parameters, such as rod placement within each bundle and rod power, which change each bundle 
40 CPR into individual rod CPR values setforth in block 707, determine the MCPR for each rod in the core, as shown 

by block 703. Individual rod power Is difficult to measure; combing that uncertainty with bundle power distribution 
uncertainty serves to increase the uncertainty in practical calculations of the SLMCPR. 

Step 4: Using the ECPR probability distribution set forth in block 708, generated by Equations 1 and 2 shown 
45 above, determine the percentage of NRSBT in the core by summing the probabilities of each rod in the core that 

is subjectto boiling transition, as shown by block 704. This summation is performed using Equation 3, shown above. 

Step 5: Vary the parameters setforth in blocks 706 and 707 for a set number of Monte Cario statistical trials, as 
shown by block 705. Compile the statistics from all the trials from steps 2 through 4 to generate a probability 
so distribution of NRSBT. 

Step 6: Compare the value of percentage of NRSBT to 0.1%, as shown in block 709. If the percentage is greater 
than 0.1%, reset the core parameters to different initial conditions in order to comply with the USNRC regulations, 
as shown in block 710. Similar to Step 1 and block 701, the new initial conditions are assumed to generate the 
55 SLMCPR. The determination of NRSBT loops until the value of NRSBT is equal to 0.1%. Similarly, if the percentage 

is less than 0.1%, the core parameters are reset to increase the value of NRSBT in order to operate the core more 
efficiently. 
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Step 7: If the percentage of NRSBT equals 0.1%, the assumed value of SLMCPR, which equals core MCPR, is 
the limit at which the core may operate, as shown by block 711 . 

Step 8: Since this process includes relatively uncertain simulations in steps 2 and 3, as shown by blocks 702 and 
5 703, the change in CPR is evaluated at a 95% confidence interval, ACPRg 5/ g 5 . The OLMCPR equals the linear 

addition of the SLMCPR to the ACPR 95/95 . The resulting value of the OLMCPR complies with the USNRC regu- 
lations. 

[0013] Motivated by the difficulties in calculating OLMCPR directly and the overly conservative approximation tech- 
10 nique currently used, the inventors were led to examine more closely some of the processes conventionally used in 
evaluating BWR designs and calculating OLMCPR. For example, the following is a brief list summarizing two of the 
most prominent factors identified by the inventors as constraining the ability to calculate OLMCPR directly using the 
ideal method: 

is 1 . The number of calculations necessary to evaluate each AOO would be too cumbersome. To establish a statis- 

tically sound estimation of the NRSBT, approximately one hundred trials for each AOO would have to be performed. 
The currently-available equipment has inherent limitations that prevent the requisite number of transient calcula- 
tions from being performed within the necessary timeframe. 

20 2. The currently-available equipment cannot simulate local power distribution (TIP uncertainty) or the relative po- 

sition of the rod within the fuel bundle (R-f actor). The variations within each rod bundle are essential to compute 
an accurate NRSBT Without a precise method of estimating the effects of the variations, a direct calculation of 
OLMCPR would be unavoidably inaccurate. 

25 [0014] Due to the excessive degree of conservatism inherent in the conventional basis used for evaluation of BWR 
operations, a substantial increase in the operating margin for a BWR can be realized by using other approaches toward 
determining the OLMCPR. Applicants have developed a mathematically-sound approach to determining OLMCPR that 
allows for a substantial increase in the permissible operating margin for a reactor and increased operating efficiency, 
greater fuel generation and/or lowered fuel exhaustion. For example, demonstrating that a particular reactor or reactor 

30 design actually has a greater operational margin than may have been previously realized, e.g., due to the use of an 
unnecessarily conservative evaluation method, could permit operation at increased output power levels or at compa- 
rable power output levels using less fuel. Accordingly, a more mathematically-sound evaluation method that results in 
a substantially greater operational margin for Boiling Water Reactors is presented herein for calculating and demon- 
strating the OLMCPR. 

35 [0015] Briefly, the improved method of the present invention is based on producing a histogram of the number of 
reactor fuel rods susceptible to operation at "boiling transition" temperatures over a range of variations in selected 
parametric quantities that are indicative of reactor design constraints and operating conditions. In addition : a core 
operational modeling approach using predetermined PDFs of the change in critical power ratio (ACPR/ICPR) for each 
type of transient event specific to the type of BWR is employed for simulating BWR thermal hydraulics and neutron 

40 kinetics during an "anticipated operational occurrence" or AOO in the reactor (for example, a operational occurrence 
that causes a brief power transient). Basically, in the present invention, all model and reactor plant parameters that 
may affect the number of rods subject to boiling transition (NRSBT) are taken into account in projecting the PDF of the 
NRSBT The NRSBT is then evaluated statistically in order to determine the OLMCPR directly without the need for first 
calculating a value for the SLMCPR. Using this approach, the present invention achieves a direct evaluation of OLMCPR 

45 of the reactor from a statistical analysis of histograms for the transient condition — as opposed the conventional "indi- 
rect" approach of determining an OLMCPR from a combination of values obtained by separately evaluating both a 
steady-state SLMCPR and the change in the critical power ratio (ACPR) due to a transient operational occurrence. 
[0016] The invention will now be described in greater detail, by way of example, with reference to the drawings, in 
which:- 

50 

Figure 1 is a schematic diagram of a simplified boiling water nuclear reactor. 

Figure 2 is a cross-sectional diagram depicting an arrangement of multiple fuel rod bundles within a core of a BWR. 
55 Figure 3 is a schematic cross-sectional diagram of an arrangement of fuel rods within a single fuel bundle. 

Figure 4 is a graph showing the determination of the NRSBT according to the ideal process. 
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Figure 5 is a flow chart illustrating a sequence of processing steps executable by a data processing system for 
performing an evaluation of OLMCPR using the ideal process. 

Figure 6 is a graph showing the linear addition of ACPR to the SLMCPR to determine OLMCPR, which is the 
currently approved process. 

Figure 7 is a flow chart illustrating a sequence of processing steps executable by a data processing system for 
performing a evaluation of the OLMCPR using the currently approved process. 

Figure 8 is a block diagram of an example data processing system used for multi-dimension/modeling of BWR 
thermal hydraulics and indirect evaluation of an OLMCPR for a BWR in accordance with the present invention. 

Figure 9 is a flow chart illustrating a sequence of processing steps used in calculating the OLMCPR using the 
generic uncertainty in ACPR/ICPR. 



Figure 1 0 is a graph showing the determination of generic uncertainty in ACPR/ICPR using the present invention. 

Figure 11 is a graph showing the determination of the NRSBT using the generic uncertainty of ACPR/ICPR. 

20 Figure 1 2 is a graph showing the comparison of the process of the present invention to the ideal process. 

Figure 13 is a graph showing the transformation of ICPR distribution to transient MCPR. 

[0017] A practical method for determining the operating limit minimum critical power ratio (OLMCPR) of a Boiling 
25 Water Reactor (BWR) nuclear core is disclosed here. This practical improvement results in the realization of greater 
operating margins for the core which in turn results in more efficient and cost effective core operation and/or configu- 
rations. This is a more direct approach towards demonstrating compliance of a nuclear reactor with USNRC licensing 
requirements than processes conventionally used for such purposes. A data processing system is disclosed including 
a computer having memory and various I/O or display devices that is specifically programmed for providing simulation 
30 of transient operational events in a BWR and for a subsequent compilation and display of one or more response 
histogram(s) that incorporate all the inherent "uncertainties" associated with reactor plant initial state conditions and 
other parameter(s) of interest or importance. 

[0018] A method is used to calculate a generic bias in change in critical power ratio during a transient event (ACPR/ 
ICPR) and uses the resulting Probability Distribution Function (PDF) to predict a more accurate OLMCPR without first 

35 calculating a SLMCPR. From a large number of experimental trials that take many factors into account, a PDF for the 
transient ACPR/ICPR is created and the standard deviation in ACPR/ICPR is determined for each transient event. A 
nominal ACPR/ICPR for the transient event starting from nominal initial conditions is also determined. Histograms of 
individual rod CPR values for the minimum point in the transient are created by drawing random values of initial CPR 
and transient ACPR/ICPR uncertainty. The initial critical power ratios (ICPR) are converted, or translated, to MCPRs 

40 by a common random value of ACPR/ICPR. From the MCPR values, the percentage of NRSBT is calculated for each 
trial. If the percentage of NRSBT is greater than 0.1%, initial operating conditions are changed and the process is 
repeated until the NRSBT is equal to 0.1%. 

[0019] The NRSBT distribution histogram is analyzed using statistical methods to determine the "central tendency" 
of the distribution. Typically the mean or median is used as a statistic to quantify central tendency. The value of this 
^5 statistic is defined here as the nominal value. In the discussions that follow, examples are given where the mean value 
is chosen as the nominal value although the present invention is not limited to this choice. Use of the median value or 
the value of some other statistic for central tendency as the nominal value is also contemplated as part of the present 
invention. 

[0020] The uncertainty in the nominal value of the statistic that is used to quantify central tendency is expressed in 
50 terms of a "confidence interval" for the nominal value. A confidence interval is defined such that there is a specified 
probability (usually of 50% or greater) that the interval contains the nominal value. For example, a 95% probability that 
the interval bounds the mean, defines a 95% confidence interval for the mean. The specified probability used to es- 
tablish this confidence interval is called the "level of confidence" or confidence level. 

[0021] The susceptibility to boiling transition during the transient is quantified statistically as either (1) the probability 
55 that a single rod in the core is susceptible to boiling transition or (2) the expected fraction of total rods in the core 
susceptible to boiling transition. Such a statistical relationship is possible because each individual trial value of NRSBT 
has been determined by summing the probabilities that individual fuel rods have CPR values less than 1 .0 during the 
transient. The nominal value for each NRSBT distribution can also by the present invention be associated with the 
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distribution of initial rod CPR values for all fuel rods in the core. It is by this process that a relationship can be established 
between the minimal initial MCPR value for all fuel rods in the core and the probability and confidence level that the 
fuel rods will be susceptible to boiling transition during the transient. The minimal initial MCPR value for the core when 
determined in this way using the probability and confidence level established by the USNRC design basis requirement 
5 for the number of rods not susceptible to boiling transition during the AOO transient, is by definition, the minimum 
Operating Limit MCPR required to demonstrate compliance. 

[0022] In accordance with one aspect, the present invention is a system including a data processing apparatus 
programmed to execute specific routines for simulating BWR core operating conditions and for calculating and statis- 
tically demonstrating the OLMCPR of a reactor in accordance with the improved method of the present invention as 

70 described in detail below. 

[0023] Figure 8 shows a block diagram of an example data processing system, contemplated for performing the 
multi-dimensional simulation of reactor core transient response and for the direct evaluation of OLMCPR for a BWR 
reactor core in accordance with the present invention. Essentially, the system includes a central processing unit 801 
(CPU), a storage memory 802, and a user interfacing I/O devices 803 and, optionally, one or more displays 804. Storage 

is memory 802 includes a data base (not shown) of reactor plant state information, parameter values and routines for 
implementing multi-dimensional simulations of core operating conditions and evaluating OLMCPR in accordance with 
the improved method of the present invention as described herein below. 

[0024] A statistical study will be performed for each type of AOO, for each class of BWR plant type, and for each 
fuel type to determine the generic transient bias and uncertainty in the ACPR/ICPR. Enough trials (on the order of one 

20 hundred) are made starting with the nominal conditions, using random variations in the model and plant parameters. 
Uncertainties in initial conditions that contribute to the ACPR/ICPR (e.g., core power) are also included in the pertur- 
bations. The data are utilized to determine bias and standard deviation on the transient ACPR/ICPR. 
[0025] A flow chart for the process of the present invention is shown in Figure 9. Block 909 remains unvaried through- 
out the calculation of the OLMCPR, and the ACPR/ICPR for individual transient events for each reactor type and fuel 

25 type must be determined before the process is used. Figure 10 shows the resulting graph of ACPR/ICPR for one 
specific type of AOO. Histogram 1000 shows the number of trials 1002 with a resulting CPR 1001 for each rod versus 
the corresponding CPR 1001 values. The PDF 1003 represents the distribution of CPR before the transient event. 
Each CPR value then changes according to individual ACPR/ICPR 1 006 values. The aggregate of the transient CPR 
values yields the PDF 1004 during the transient event The nominal ACPR/ICPR 1005 is defined to be the difference 

30 in the nominal CPR value of the PDF 1 003 and the nominal CPR value of the PDF 1 004. The calculation of the OLMCPR 
is as follows. 



Step 1 : Assume a set of base core operating conditions using the parameters to run the plant generates a core 
MCPR equal to the OLMCPR as shown by block 901 . 



35 



Step 2: Using the parameters, such as core power, core flow, core pressure, bundle power and others, that predict 
a general bundle CPR set forth in block 907, determine the ICPR for each bundle in the core, as shown by block 902. 



40 



Step 3: Using parameters, such as rod placement within each bundle and rod power distribution, that change each 
bundle CPR into individual rod CPR values set forth in block 908, determine the ICPR for each rod in the core, as 
shown by block 903, 



45 



Step 4: Using a randomly drawn individual ACPR/ICPR 1006 value from the graph of the appropriate transient 
represented in Figure 10, MCPR values are projected for corresponding values of ICPR according to Equation 4. 
In Figure 11, this process is represented by Shift 1109. Histogram 1100 shows the number of rods at a specific 
CPR value 1102 versus the corresponding CPR value 1101. The histogram 1107 is translated to histogram 1108 
during the transient using a randomly selected ACPR/ICPR 1 006 value. Lowest CPR value 1 1 05 becomes lowest 
CPR value 1106, and lowest CPR rod PDF 1103 becomes lowest CPR rod 1104. 



50 




55 



(Equation 4) 



Step 5: Using the ECPR probability distribution shown as PDF 1104 and set forth in block 910, determine the 



7 



BNSDOCID: <EP 1221701A1_I_> 



10 



15 



25 



30 



EP 1 221 701 A1 

percentage of NRSBT in the core by summing the probabilities of each rod in the core that is subject to boiling 
transition as shown by block 905. This summation is performed using Equation 3, shown above. 

Step 6: Vary the parameters set forth in blocks 907 and 908 for a set number of Monte Carlo statistical trials as 
shown by block 906. Compile the statistics from all the trials from steps 2 through 5 to qenerate a probability 
distribution of NRSBT 

Step 7: Compare the value of percentage of NRSBT to 0.1% as shown in block 911. If the percentage is greater 
than 0.1%, reset the core parameters to different initial conditions in order to comply with the USNRC regulations 
as shown in block 912. Similar to Step 1 and block 901 , the new initial conditions are assumed to generate the 
OLMCPR. The determination of NRSBT restarts and runs until the value of NRSBT is equal to 0.1%. Similarly if 
the percentage is less than 0.1%, the core parameters are reset to increase the value of NRSBT in order to operate 
the core more efficiently or to reduce effluents. 

Step 8: If the percentage of NRSBT equals 0.1%, the assumed value of OLMCPR, which equals core MCPR 
complies with the USNRC regulations as shown by block 913. Accordingly, the operating core conditions are set 
as the assumed parameters. 

[0026] Two assumptions are made for the above estimation of OLMCPR. First, in performing step 4, shown in Figure 
1 1 as shift 1 1 09 and in Figure 9 as block 904, the inventors assume that random draws from the ACPR/ICPR distribution 
are permissible for a perturbation in the initial conditions. Therefore, variations in ACPR/ICPR must be independent of 
perturbations in initial conditions or have a negative correlation, so that the interaction tends to diminish the individual 
effects. Second, in performing step 4, the inventors assume that the transient change in the ACPR/ICPR applies to all 
rods. rr 

[0027] A demonstration analysis shows that the ACPR/ICPR is not sensitive to the uncertainty in core power core 
flow, core pressure, feedwater temperature, and rod peaking factor (R-factor). Of these, one of the most important 
parameters in the currently approved process is core power. This parameter actually results in an effect opposite the 
effect on ICPR. If the power increases, the ICPR will decrease but the ACPR/ICPR will also decrease This will result 
in an MCPR that would be higher than derived through the currently approved process. Another conservative factor is 
the intended use of the nominal ACPR/ICPR. If the core was adjusting to a limiting rod pattern to maximize the number 
of contributing bundles, as is done for the currently approved process, the ACPR/ICPR is 4% lower. 

IMPACT OF CRITICAL ICPR UNCERTAINTIES ON ACPR/ICPR (Table 1) 



35 


Quantity 


Uncertainty in 
Quantity 


ACPR/ICPR impact at 1a (%) 




Total Core Power 


+1a = +2% 
-1a = -2% 


-0.6% 
+0.7% 


40 


Total Core Flow 


+1a = +2.5% 
-1a = -2.5% 


+4.4% 
-0.2% 




Steam Dome Pressure 


+1a = +0.7% 
-1a = -0.7% 


-0.5% 
+0.6% 


45 


Feedwater Temperature 


+1a = +0.8% 
-1a = -0.8% 


+0.3% 
-0.3% 




R-factor 


+1a = +1.6% 
-1a- -1.6% 


-0.5% 
+0.5% 


50 


Radial Power Distribution (increase radial peaking with 
SLMCPR type limiting rod pattern) 


+1a = 3.19% 


-5.0% 




(column 101) 


(column 102) 


(column 103) 



55 



[0028] Table 1 shows the impact of uncertainty in critical ICPR values on ACPR/ICPR values. Column 101 lists the 
critical parameter quantities that affect the ACPR/ICPR. Column 102 lists the percentage uncertainty of each parameter 
corresponding to the standard deviation of the associated PDF. a is the standard deviation of the PDF corresponding 
to the uncertainty in parameter quantity. Column 103 lists the change in the ACPR/ICPR corresponding to a change 
of one standard deviation of each parameter. 
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[0029] The ACPR/ICPR is not sensitive to the other unknown parameter in the currently approved process. The axial 
power distribution is also part of the local power distribution (TIP uncertainty) calculation in the currently approved 
process. For a very large change in axial power shape (nearly two times higher power in the bottom of the bundle), 
the sensitivity to ACPR/ICPR is less than 2%, which is insignificant relative to TIP uncertainty. 
5 [0030] The other assumption to be validated is that a constant value of ACPR/ICPR can be applied to rods at different 
ICPR values. As described above, the transient MCPR distribution will be obtained by transforming the ICPR distribution 
using Equation 4. 

[0031] To further validate this assumption, a specific set of calculations were performed. Benchmark calculations 
were made for a transient event that included the uncertainties in core power and channel pressure drop as initial 

10 conditions, as well as uncertainties in the model. Core power and channel pressure drop uncertainties were chosen, 
because they are the only currently approved process compatible uncertainties that are also varied in generating the 
generic uncertainty probability distribution function. MCPR distributions during the transient were generated for two 
fuel bundles in the core through ninety-eight transient calculations. The two bundles are very close in ICPR values and 
. have identical ACPR/ICPR values. 

15 [0032] To verify the translation process, ninety-eight Monte Carlo calculations were then performed where only the 
core power and pressure drop were varied to generate a PDF of ICPRs at the initial operating state. Figure 12 shows 
histogram 1200, which is the number of rods 1202 at a certain CPR versus the corresponding CPR 1201 value. PDF 
1203 is the ICPR distribution that was created using the Monte Carlo calculations varying core power and pressure 
drop. PDF 1205 is the corresponding transient MCPR distribution after the process of the invention transformation was 

20 applied. PDF 1 204 is the reference ICPR distribution. PDF 1 206 is the transient MCPR distribution when applying the 
currently approved process. PDF 1205 and PDF 1206 are very similar in both the most probably value of MCPR and 
the associated standard deviation of each distribution. Since there is a strong resemblance between the two resulting 
MCPR distributions, the transformation using the process of the invention is valid. 

[0033] It has been demonstrated that: (1) the ACPR/ICPR is independent relative to the uncertainties that affect the 
25 ICPR, or the covariance is such that it is conservative to assume independence and (2) the transient MCPR distribution 
can be determined by applying the transient ACPR/ICPR uncertainty to the rod ICPR distribution using the proposed 
approach. 

[0034] An example of the process of the invention is described by Figure 13. In Figure 13, histogram 1300 shows 
the number of rods 1302 of a certain CPR value versus the corresponding CPR value 1301 . The PDF 1303 shows the 

30 resulting ICPR values from a set of approximately ninety-eight ICPR trials with all uncertainties applied. Ninety-eight 
new trials were run to generate a ACPR/ICPR distribution for the specific transient event in order to translate the ICPR 
values to MCPR values. This ACPR/ICPR distribution is not shown in Figure 13. The ACPR/ICPR distribution was 
applied using the process of the invention to the ICPR PDF 1303 to obtain the MCPR PDF 1304. The NRSBT was 
determined using the process of the invention, and the OLMCPR was determined to be 1 .26. As a comparison, using 

35 the currently approved process, the SLMCPR was determined to be 1 .1 0. Thus, the process described herein is more 
conservative than the first stage of the currently approved process. Ultimately, however, the currently approved process 
generates a unnecessarily conservative value after the error factor is added to the SLMCPR value, which yields a 
OLMCPR value needlessly larger than the process of the invention. 

[0035] Although the improved methods, as described herein below, are preferably implemented using a high speed 
^o data processing system capable of processing simulation routines that require highly accurate calculations and multiple 
reiterations, the present invention is not intended as limited to any one particular type of computer or data processing 
system. Any generic data processing system having sufficient, speed, storage memory and programmable computa- 
tional capabilities for implementing statistical data analysis/reduction may be utilized to implement the present inven- 
tion. 

[0036] For the sake of good order, various aspects of the invention are set out in the following clauses:- 

1. A method for direct evaluation of an operating limit minimum critical power ratio (OLMCPR) of a boiling water 
reactor (BWR) using computer simulation (s) of postulated operational events and previously determined transla- 
tions in change in critical power ratio (ACPR/ICPR) for transient events, said reactor described by one or more 
50 parametric quantities indicative of design constraints and operating conditions, comprising the steps: 

a) determining a probability distribution P(ACPR/ICPR) of potential ACPR/ICPR for the class of BRW, opera- 
tional event, and fuel type; 

b) initializing said parametric quantities to nominal values; 

55 c) determining steady state initial conditions of reactor operation; 

d) simulating an operational event for a plurality of fuel rods in a BWR core; 

e) calculating a minimum critical power ratio (MCPR) for each fuel rod simulated in step (d) using the probability 
distribution obtained in step (a); 
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f) determining a probability distribution P(MCPR) of potential MCPR values for each MCPR value obtained in 
step (e); 

g) calculating a value for a total number of fuel rods subject to boiling transition (NRSBT) by summing portions 
of probability distributions corresponding to values for MCPR where MCPR<1 .0 for each probability distribution 
from step (f); 

h) perturbing one or more of said parametric quantities and recalculating another value for NRSBT; 

i) repeating steps (c) through (h) for a predetermined number of perturbations; 
j) developing a histogram of NRSBT values calculated in steps (g) through (h); 

k) calculating a nominal NRSBT value, based on a central tendency of an NRSBT distribution, from the his- 
togram of NRSBT values compiled in step (j); 
I) calculating a confidence interval for the nominal NRSBT value; 

m) selecting an OLMCPR for the reactor as the initial minimal MCPR such that, for a simulation of the most 
limiting operational event the nominal NRSBT value at a prescribed level of confidence remains less than a 
predetermined cutoff value; and 

n) effecting said BWR operation by applying the OLMCPR selected in step (I) as an operational control pa- 
rameter. 



35 



2. The method of clause 1 wherein the operational event is a transient event associated with an Anticipated Op- 
erational Occurrence (AOO). 

20 

3. The method of clause 1 , wherein multi-dimensional modeling of reactor thermal hydraulics and power is used 
to simulate an anticipated operational occurrence. 

4. The method of clause 1 , wherein the parametric quantities correspond to reactor plant state values and/or 
25 modeling parameters. 

5. The method of clause 1 , wherein step (d) the anticipated occurrence is simulated for a plurality of fuel rods 
simultaneously. 

30 6. The method of clause 1 , wherein said perturbing of one or more of said parametric quantities is accomplished 

using a Monte-Carlo statistical analysis approach. 

7. A system for determining an operating limit minimum critical power ratio (OLMCPR) of a boiling water reactor 
(BWR), said system comprising a computer including a storage memory and I/O devices, said memory having 
stored therein rules for simultaneously simulating and evaluating thermal operating characteristics for a plurality 
of fuel rods during an anticipated operational occurrence of the reactor and a data base of one or more parametric 
quantities representing reactor plant operational state values and/or fuel rod modeling parameters, said computer 
programmed to: 

40 (') simulate an anticipated operational occurrence for the reactor; 

(ii) determine a minimum critical power ratio (MCPR) for each individual fuel rod simulated based on a prede- 
termined probability distribution P(ACPR/ICPR) of potential change in critical power ratio (ACPR/ICPR) values 
for said anticipated operational occurrence; 

(iii) calculate a probability distribution P(MCPR) of potential MCPR values for each MCPR value; 

(iv) calculate a value for a number of fuel rods subject to boiling transition (NRSBT) by summing portions of 
probability distributions corresponding to values for MCPR where MCPR<1.0 for all rods simulated, then re- 
petitively calculate further values for NRSBT after perturbing one or more of said parametric quantities for a 
predetermined number of different perturbations; 

(v) develop a histogram of calculated NRSBT values over all perturbations; and 

(vi) select the OLMCPR from a minimal MCPR such that, for a simulation of the most limiting operational event, 
the nominal NRSBT value at a prescribed level of confidence as determined from analysis of the histogram 
determined from step (v) remains less than a predetermined cutoff value. 

8. The system of clause 7 wherein the operational event simulated is a transient event associated with an Antici- 
5 $ pated Operational Occurrences (AOO) 

9. The system of clause 7 wherein a multi-dimensional modeling of fuel rod thermal hydraulics and reactor power 
is used to simulate an anticipated operational occurrence. 



45 



50 
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10. The system of clause 7 wherein step said anticipated occurrence is simulated for a plurality of fuel rods simul- 
taneously. 

11. A computer program product embodied on a computer readable medium for determining an operating limit 
5 minimum critical power ratio (OLMCPR) for a Boiling Water Reactor (BWR), comprising: 

first instruction sequence means for performing a multi-dimensional simulation of reactor thermal hydraulics 
and power during an operational event in the reactor using a predetermined probability distribution of change 
in critical power ratio (ACPR/ICPR) for said operational event; and 
10 second instruction sequence means, coupled to said first instruction means, for directly calculating an OLM- 

CPR value based on a statistic evaluation of a minimum critical power ratio (MCPR) for each fuel rod, as 
obtained from said first instruction sequence means, and for developing a histogram of possible NRSBT 
(number of fuel rods subject to boiling transition) values calculated from perturbations of predetermined reactor 
plant state values and fuel rod modeling parameters. 

75 

12. The computer program of clause 11 wherein said means for performing a multi-dimensional simulation of 
reactor thermal hydraulics and power simulates a transient operational occurrence for a plurality of fuel rods si- 
multaneously. 

20 13. A method for evaluating the operating limit minimum critical power ratio (OLMCPR) of a boiling water reactor 

(BWR), comprising the steps, executed by a computer, of: 

a) setting a minimum critical power ratio (MCPR) for each of a plurality of individual nuclear fuel rods used in 
the reactor; 

25 b) calculating a value for a total number of fuel rods subject to boiling transition (NRSBT) based on a pre- 

determined distribution of change in critical power ratio (ACPR/ICPR); 

c) perturbing one or more reactor plant state value and/or fuel rod modeling parameter and recalculating a 
value for NRSBT; 

d) performing step (c) for a predetermined number of perturbations; 

30 e) developing a histogram of NRSBT values determined in steps (b), (c) and (d); and 

f) selecting the OLMCPR from a minimal MCPR such that, for a simulation of the most limiting operational 
event, the nominal NRSBT value at a prescribed level of confidence as determined from analysis of the his- 
togram determined from step (e) remains less than a predetermined cutoff value. 

35 14. The method of clause 1 3 wherein the calculating of an NRSBT value in step (b) is accomplished by the steps of: 

determining a probability distribution P(MCPR) to each MCPR value set in step (a), said P(MCPR) indicative 
of a range of possible MCPR values resulting from various operational and design uncertainties; and 
integrating the probability distribution for values of MCPR<1 .0 for each fuel rod and summing the integration 
40 results for all fuel rods. 

15. The method of clause 13 wherein the perturbing in step (c) of one or more reactor plant state value(s) and/or 
modeling parameter(s) is accomplished using randomly generated variations of said values and parameters. 

45 16. For a nuclear fuel core of a Boiling Water Reactor, wherein fuel design and/or core configuration are contingent 

upon an operating margin for the reactor, said operating margin being determined by a process for evaluating an 
operating limit minimum critical power ratio (OLMCPR), an improved process for evaluating an OLMCPR and 
operating a Boiling Water Reactor which results in an improved operating margin for the reactor, comprising the 
steps of: 

50 

a) developing in a memory of a computer system a histogram of N RSBT (number of fuel rods subject to boiling 
transition) values, corresponding to a plurality of computer simulations of a transient operational occurrence 
in a reactor based on a predetermined probability distribution of the change in critical power ratio (ACPR/ 
ICPR), said simulations providing values for a critical power ratio (CPR) for one or more fuel rods for a multiple 

55 of different parametric quantities for said reactor; 

b) selecting a nominal NRSBT value, based on a central tendency of NRSBT distribution, statistically deter- 
mined by said computer system from the histogram of NRSBT values developed in step (a); 

c) selecting a confidence interval for the nominal NRSBT value; 
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d) selecting an OLMCPR value from a minimal CPR such that during a simulation of a limiting transient oper- 
ational occurrence the nominal NRSBT value remains less than a predetermined cutoff value; and 

e) effecting said Boiling Water Reactor operation by applying the OLMCPR selected in step (d) as an opera- 
tional control parameter. 

17. The method clause 16 wherein said histogram is developed on a data processing system using multi-dimen- 
sional modeling of transient operational occurrences, said data processing system including memory to store CPR 
data obtained from said simulations. 

18. For a nuclear fuel core of a Boiling Water Reactor, wherein fuel rod design and/or core configuration are con- 
tingent upon an operating margin for the reactor, said operating margin being determined by a process for evalu- 
ating an operating limit minimum critical power ratio (OLMCPR), a improved process for evaluating an OLMCPR 
and operating a Boiling Water Reactor that results in an increased operating margin for the reactor, comprising 
the steps, executed by a computer, of: 

a) developing in a memory of said computer a histogram of NRSBT (number of fuel rods subject to boiling 
transition) values, corresponding to a plurality of computer simulations of a transient operational occurrence 
in a reactor, said simulations providing values for a critical power ratio (CPR) for one or more fuel rods for a 
multiple of different parametric quantities for said reactor, said simulations based on a predetermined proba- 
bility distribution of change in critical power ratio (ACPR/ICPR) for said transient operational occurrence; 

b) calculating a nominal NRSBT value, based on a central tendency of an NRSBT distribution, from the his- 
togram of NRSBT values obtained in step (a): 

c) selecting a confidence interval for the nominal NRSBT value; 

d) selecting an OLMC PR for the reactor as the initial minimal MCPR such that during a simulation of a transient 
the nominal NRSBT value remains less than a predetermined cutoff value; and 

e) effecting Boiling Water Reactor operation by applying the OLMCPR selected in step (d) as an operational 
control parameter. 

19. A method for statistically demonstrating an operating limit minimum critical power ratio (OLMCPR) of a boiling 
water reactor (BWR) for compliance with licensing requirements, said reactor characterized by one or more par- 
ametric quantities indicative of design constraints and operating conditions, comprising the steps of: 

a) programming a computer to determine an OLMCPR value for a boiling water reactor, said computer pro- 
grammed at least to: 

i) develop a histogram of NRSBT (number of fuel rods subject to boiling transition) values, corresponding 
to a plurality of simulations of a transient operational occurrence in a reactor, said simulations providing 
values for a critical power ratio (CPR) for one or more fuel rods over a multiple of randomly selected 
variations in the parametric quantities for said reactor, said simulations based on a predetermined prob- 
ability distribution of change in critical power ratio (ACPR/ICPR) for said transient operational occurrence; 

ii) calculate a nominal NRSBT value from the histogram of NRSBT values; and 

iii) select an OLMCPR for the reactor corresponding to a minimal MCPR such that during a simulation of 
a transient the nominal NRSBT value remains less than a predetermined cutoff value; and 

b) providing said OLMCPR value to an output device for display recordation or storage. 

20. For use in controlling a Boiling Water Reactor having a nuclear fueled core characterized by an operating limit 
minimum critical power ratio (OLMCPR) value, a method for determining said OLMCPR value comprising the steps 
of: 

a) using a computer to simulate transient operational occurrences based on predetermined distributions of 
said occurrences which might occur during the operation of a Boiling Water Reactor; 

b) developing a histogram of NRSBT values from a computer simulation of a transient operational occurrence 
in a Boiling Water Reactor, the simulation providing values for a critical power ratio (CPR) for one or more fuel 
rods; 

c) determining a nominal NRSBT value, based on a central tendency of an NRSBT distribution, from the 
histogram of NRSBT values obtained in step (a); 

d) selecting a confidence interval for the nominal NRSBT value; and 
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e) selecting an OLMCPR for the reactor as the initial minimal MCPR such that during a simulation of a transient 
the nominal NRSBT value remains less than a predetermined cutoff value. 



s Claims 

1. A method for direct evaluation of an operating limit minimum critical power ratio (OLMCPR) of a boiling water 
reactor (BWR) using computer simulation(s) of postulated operational events and previously determined transla- 
tions in change in critical power ratio (ACPR/ICPR) for transient events, said reactor described by one or more 

10 parametric quantities indicative of design constraints and operating conditions, comprising the steps: 

a) determining a probability distribution P(ACPR/ICPR) of potential ACPR/ICPR for the class of BRW, opera- 
tional event and fuel type; 

b) initializing said parametric quantities to nominal values; 

15 c) determining steady state initial conditions of reactor operation; 

d) simulating an operational event for a plurality of fuel rods in a BWR core; 

e) calculating a minimum critical power ratio (MCPR) for each fuel rod simulated in step (d) using the probability 
distribution obtained in step (a); 

f) determining a probability distribution P(MCPR) of potential MCPR values for each MCPR value obtained in 
20 step(e); 

g) calculating a value for a total number of fuel rods subject to boiling transition (NRSBT) by summing portions 
of probability distributions corresponding to values for MC PR where MCPR<1 .0 for each probability distribution 
from step (f); 

h) perturbing one or more of said parametric quantities and recalculating another value for NRSBT; 
25 i) repeating steps (c) through (h) for a predetermined number of perturbations; 

j) developing a histogram of NRSBT values calculated in steps (g) through (h); 

k) calculating a nominal NRSBT value, based on a central tendency of an NRSBT distribution, from the his- 
togram of NRSBT values compiled in step (j); 
I) calculating a confidence interval for the nominal NRSBT value; 
30 m) selecting an OLMCPR for the reactor as the initial minimal MCPR such that, for a simulation of the most 

limiting operational event the nominal NRSBT value at a prescribed level of confidence remains less than a 
predetermined cutoff value; and 

n) effecting said BWR operation by applying the OLMCPR selected in step (I) as an operational control pa- 
rameter. 

35 

2. The method of claim 1 wherein the operational event is a transient event associated with an Anticipated Operational 
Occurrence (AOO). 

3. A system for determining an operating limit minimum critical power ratio (OLMCPR) of a boiling water reactor 
40 (BWR), said system comprising a computer including a storage memory and I/O devices, said memory having 

stored therein rules for simultaneously simulating and evaluating thermal operating characteristics for a plurality 
of fuel rods during an anticipated operational occurrence of the reactor and a data base of one or more parametric 
quantities representing reactor plant operational state values and/or fuel rod modeling parameters, said computer 
programmed to: 

45 

(i) simulate an anticipated operational occurrence for the reactor; 

(ii) determine a minimum critical power ratio (MCPR) for each individual fuel rod simulated based on a prede- 
termined probability distribution P(ACPR/ICPR) of potential change in critical power ratio (ACPR/ICPR) values 
for said anticipated operational occurrence; 

50 (iii) calculate a probability distribution P(MCPR) of potential MCPR values for each MCPR value; 

(iv) calculate a value for a number of fuel rods subject to boiling transition (NRSBT) by summing portions of 
probability distributions corresponding to values for MCPR where MCPR<1 .0 for all rods simulated, then re- 
petitively calculate further values for NRSBT after perturbing one or more of said parametric quantities for a 
predetermined number of different perturbations; 

55 (v) develop a histogram of calculated NRSBT values over all perturbations; and 

(vi) select the OLMCPR from a minimal MCPR such that, for a simulation of the most limiting operational event, 
the nominal NRSBT value at a prescribed level of confidence as determined from analysis of the histogram 
determined from step (v) remains less than a predetermined cutoff value. 
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4. The system of claim 3 wherein the operational event simulated is a transient event associated with an Anticipated 
Operational Occurrences (AOO) 

5. A computer program product embodied on a computer readable medium for determining an operating limit mini- 
mum critical power ratio (OLMCPR) for a Boiling Water Reactor (BWR), comprising: 

first instruction sequence means for performing a multi-dimensional simulation of reactor thermal hydraulics 
and power during an operational event in the reactor using a predetermined probability distribution of change 
in critical power ratio (ACPR/ICPR) for said operational event; and 

second instruction sequence means, coupled to said first instruction means, for directly calculating an OLM- 
CPR value based on a statistic evaluation of a minimum critical power ratio (MCPR) for each fuel rod, as 
obtained from said first instruction sequence means, and for developing a histogram of possible NRSBT 
(number of fuel rods subject to boiling transition) values calculated from perturbations of predetermined reactor 
plant state values and fuel rod modeling parameters. 

6. A method for evaluating the operating limit minimum critical power ratio (OLMCPR) of a boiling water reactor 
(BWR), comprising the steps, executed by a computer, of: 



a) setting a minimum critical power ratio (MCPR) for each of a plurality of individual nuclear fuel rods used in 
20 the reactor; 

b) calculating a value for a total number of fuel rods subject to boiling transition (NRSBT) based on a pre- 
determined distribution of change in critical power ratio (ACPR/ICPR); 

c) perturbing one or more reactor plant state value and/or fuel rod modeling parameter and recalculating a 
value for NRSBT; 

25 d) performing step (c) for a predetermined number of perturbations; 

e) developing a histogram of NRSBT values determined in steps (b), (c) and (d); and 

f) selecting the OLMCPR from a minimal MCPR such that, for a simulation of the most limiting operational 
event, the nominal NRSBT value at a prescribed level of confidence as determined from analysis of the his- 
togram determined from step (e) remains less than a predetermined cutoff value 

30 

7. For a nuclear fuel core of a Boiling Water Reactor, wherein fuel design and/or core configuration are contingent 
upon an operating margin for the reactor, said operating margin being determined by a process for evaluating an 
operating limit minimum critical power ratio (OLMCPR), an improved process for evaluating an OLMCPR and 
operating a Boiling Water Reactor which results in an improved operating margin for the reactor, comprising the 
35 steps of: 

a) developing in a memory of a computer system a histogram of NRSBT (number of fuel rods subject to boiling 
transition) values, corresponding to a plurality of computer simulations of a transient operational occurrence 
in a reactor based on a predetermined probability distribution of the change in critical power ratio (ACPR/ 

40 ICPR), said simulations providing values for a critical power ratio (CPR) for one or more fuel rods for a multiple 

of different parametric quantities for said reactor; 

b) selecting a nominal NRSBT value, based on a central tendency of NRSBT distribution, statistically deter- 
mined by said computer system from the histogram of NRSBT values developed in step (a); 

c) selecting a confidence interval for the nominal NRSBT value; 

45 d ) selecting an OLMCPR value from a minimal CPR such that during a simulation of a limiting transient oper- 

ational occurrence the nominal NRSBT value remains less than a predetermined cutoff value; and 
e) effecting said Boiling Water Reactor operation by applying the OLMCPR selected in step (d) as an opera- 
tional control parameter. 



8. A process for evaluating an operating limit minimum critical power ratio (OLMCPR) for a nuclear fuel core of a 
Boiling Water Reactor that results in an increased operating margin for the reactor, wherein fuel rod design and/ 
or core configuration are contingent upon an operating margin for the reactor, said operating margin being deter- 
mined by a process for evaluating OLMCPR, the process comprising the steps, executed by a computer, of: 

a) developing in a memory of said computer a histogram of NRSBT (number of fuel rods subject to boiling 
transition) values, corresponding to a plurality of computer simulations of a transient operational occurrence 
in a reactor, said simulations providing values for a critical power ratio (CPR) for one or more fuel rods for a 
multiple of different parametric quantities for said reactor, said simulations based on a predetermined proba- 
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bility distribution of change in critical power ratio (ACPR/ICPR) for said transient operational occurrence; 

b) calculating a nominal NRSBT value, based on a central tendency of an NRSBT distribution, from the his- 
togram of NRSBT values obtained in step (a); 

c) selecting a confidence interval for the nominal NRSBT value; 

5 d) selecting an OLMCPR for the reactor as the initial minimal MCPR such that during a simulation of a transient 

the nominal NRSBT value remains less than a predetermined cutoff value; and 

e) effecting Boiling Water Reactor operation by applying the OLMCPR selected in step (d) as an operational 
control parameter. 

10 9. A method for statistically demonstrating an operating limit minimum critical power ratio (OLMCPR) of a boiling 
water reactor (BWR) for compliance with licensing requirements, said reactor characterized by one or more par- 
ametric quantities indicative of design constraints and operating conditions, comprising the steps of: 

a) programming a computer to determine an OLMCPR value for a boiling water reactor, said computer pro- 
's grammed at least to: 

i) develop a histogram of NRSBT (number of fuel rods subject to boiling transition) values, corresponding 
to a plurality of simulations of a transient operational occurrence in a reactor, said simulations providing 
values for a critical power ratio (CPR) for one or more fuel rods over a multiple of randomly selected 

20 variations in the parametric quantities for said reactor, said simulations based on a predetermined prob- 

ability distribution of change in critical power ratio (ACPR/ICPR) for said transient operational occurrence; 

ii) calculate a nominal NRSBT value from the histogram of NRSBT values; and 

Hi) select an OLMCPR for the reactor corresponding to a minimal MCPR such that during a simulation of 
a transient the nominal NRSBT value remains less than a predetermined cutoff value; and 



25 



b) providing said OLMCPR value to an output device for display, recordation or storage. 



10. For use in controlling a Boiling Water Reactor having a nuclear fueled core characterized by an operating limit 
minimum critical power ratio (OLMCPR) value, a method for determining said OLMCPR value comprising the steps 
30 of: 

a) using a computer to simulate transient operational occurrences based on predetermined distributions of 
said occurrences which might occur during the operation of a Boiling Water Reactor; 

b) developing a histogram of NRSBT values from a computer simulation of a transient operational occurrence 
35 jn a Boiling Water Reactor, the simulation providing values for a critical power ratio (CPR) for one or more fuel 

rods; 

c) determining a nominal NRSBT value, based on a central tendency of an NRSBT distribution, from the 
histogram of NRSBT values obtained in step (a); 

d) selecting a confidence interval for the nominal NRSBT value; and 

40 e) selecting an OLMC PR for the reactor as the initial minimal MCPR such that during a simulation of a transient 

the nominal NRSBT value remains less than a predetermined cutoff value. 
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